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Abstract. 

\mJ ' An anomalously large transparency of the Universe to gamma rays has recently been discovered 

by the Imaging Atmospheric Cherenkov Telescopes (IACTs) H.E.S.S. and MAGIC. We show that 
observations can be reconciled with standard blazar emission models provided photon oscillations 
into a very light Axion-Like Particle occur in extragalactic magnetic fields. A quantitative estimate 
of this effect is successfully applied to the blazar 3C279. Our prediction can be tested with the 
satellite-borne FermifLNI detector as well as with the ground-based IACTs H.E.S.S., MAGIC, 
CANGAROO ITT, VERITAS and the Extensive Air Shower arrays ARGO-YBT and MILAGRO. Our 
result also offers an important observational test for models of dark energy wherein quintessence is 
coupled to the photon through an effective dimension-five operatofl 
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An impressive amount of information about the Universe in the very-high-energy (VHE) 
band has been collected over the last few years by the Imaging Atmospheric Cherenkov 
Telescopes (IACTs) H.E.S.S., MAGIC, CANGAROO III and VERITAS. Observations 
carried out by these IACTs concern gamma-ray sources over an extremely wide interval 
of distances, ranging from the parsec scale for Galactic objects up to the Gigaparsec 
scale for the farthest detected blazar 3C279. This circumstance allows not only to infer 
the intrinsic properties of the sources, but also to probe the nature of photon propagation 
throughout cosmological distances. 

The latter isssue becomes particularly important in the energy range above lOOGeV, 
where the horizon of the observable Universe rapidly shrinks as the energy further in- 
creases. This comes about because photons from distant sources scatter off background 
photons permeating the Universe, thereby disappearing into electron-positron pairs |Q]]. 
The corresponding cross section d(yy — > e + e~) turns out to peak when the VHE photon 
energy E and the background photon energy e are related by £ ~ (500GeV /E) eV, 
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so that the resulting cosmic opacity is dominated by the interaction with ultravio- 
let/optical/infrared photons of the diffuse extragalactic background - usually called ex- 
tragalactic background light (EBL) - for observations performed by IACTs. 

Owing to the absorption process in question, photon propagation is controlled by the 
optical depth t(E,D), with D denoting the source distance. Therefore, the observed 
photon spectrum 4> t, s (.E,D) is related to the emitted one <& em (E) by 

^ ohs (E,D)=e-^ E ' D ^ em (E). (1) 

Unlike the CMB, the EBL is produced by galaxies during the whole age of the Universe 
and possibly by a first generation of stars formed before galaxies were assembled. 
Based on stellar evolution models in galaxies as well as on deep galaxy counts, several 
groups have attempted a determination of the spectral energy distribution of the EBL 
and ultimately of the optical depth x{E,D) for 100 GeV < E < lOOTeV [@]. Because 
galaxies were brigther in the past than they are now, evolutionary effects should in 
principle be included in the evaluation of t(E,D) but they become unimportant at 
sufficiently low redshift and will therefore be neglected throughout. Correspondingly 
we have t(E,D) ~ D/Xy(E), with Xy(E) denoting the photon mean free path for yy — > 
e + e~. As a consequence, Eq. (Q} becomes 

4> obs (£,D) ^ e~ D M E ) $ em (£) . (2) 

The function Xy(E) has been computed within realistic models of the EBL and is 
reported in Fig. 1 (from ref. [0]). 

We see that Xy{E) decreases roughly like a power law from the Hubble radius 4.3 Gpc 
slightly below lOOGeV to about IMpc at lOOTeV. Thus, Eq. © entails that the ob- 
served flux is exponentially suppressed both at high energy and at large distances, so 
that sufficiently far-away sources become hardly visible in the VHE range. Moreover, 
the observed spectrum of distant sources gets much steeper than the emitted one. 

Yet, observations have not detected the behaviour predicted by Eq. ©. A first indi- 
cation in this direction was reported by the H.E.S.S. collaboration in connection with 
the discovery of the two blazars H2356-309 (z = 0.165) and 1ES 1101-232 (z = 0.186) 
at E ~ ITeV [4]. Stronger evidence comes from the observation of the blazar 3C279 
(z = 0.536) at E ~ 0.5 TeV by the MAGIC collaboration [5]. In particular, the signal 
from 3C279 collected by MAGIC in the region E < 220 GeV has more or less the same 
statistical significance as the one in the range 220 GeV < E < 600 GeV (6.1c in the 
former case, 5. 1 a in the latter) dH]. 

Turning the argument around and assuming standard photon propagation, the ob- 
served spectrum can only be reproduced by an emission spectrum much harder than for 
any other observed blazar. Such a hard spectrum is also largely inconsistent with the 
predictions of current blazar models @]. So, we are led to believe that the effect under 
consideration ought to be due to an anomalous photon propagation over cosmological 
distances rather than to new phenomena occurring inside the blazars themselves. From 
this point of view, the Universe appears to be more transparent to VHE gamma rays 
than previously thought, up to the point that a conflict with the standard scenario can be 
envisaged. 
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FIGURE 1. Pair-production photon mean free path as a function of beam photon energy (fromref. H). 



A cautionary remark is in order. Available observational information is insufficient to 
fully settle the issue and in fact it has been pointed out that an anomalously hard emission 
spectrum can be produced within unconventional blazar models [7] . However, according 
to this option it looks difficult to explain why the physical mechanisms occurring in the 
most distant blazars are so different from those described by standard blazar models. 

A way out of the considered difficulty has recently been proposed by the present 
authors and rests upon an oscillation mechanism occurring in extragalactic magnetic 
fields, whereby a photon can become a new very light spin-zero particle, named Axion- 
Like Particle (ALP) [8]. Once produced, the ALP travels unimpeded throughout the 
Universe and can convert back to a photon before reaching the Earth, thereby acting 
as if the observed photons had travelled a distance largely in excess of their mean free 
path. Thanks to the exponential dependence of the observed flux on Xy(E), even a small 
increase in Xy(E) gives rise to a large effect. More specifically, this mechanism yields 
an observed spectrum much harder than the one predicted by Eq. © for distant blazars, 
thereby leading to agreement with observations for standard emission spectra. 

Our aim is to review the main features of our proposal as well as its application to 
blazar 3C 279. 



PHOTON-ALP OSCILLATIONS 



Motivations for ALPs 

In spite of the enormous success scored by the Standard Model of strong, weak and 
electromagnetic interactions in describing physics at energies up to the Fermi scale 

— 1/2 

G F ~ 250 GeV, nobody would seriously regard it as the final theory. Instead, on 
the basis of phenomenological as well as conceptual reasons, the Standard Model is 
currently viewed as the low-energy manifestation of some more fundamental and richer 
theory of all elementary-particle interactions including gravity. Therefore, the lagrangian 
of the Standard Model is expected to be modified by small non-renormalizable terms 
describing interactions among known and new particles. 

Attempts to go beyond the Standard Model are the central research topic in high- 
energy physics since almost three decades and each specific proposal is characterized 
by a set of new particles along with their mass spectrum and interactions with the 
standard world. ALPs turn out to be a generic prediction of many extensions of the 
Standard Model and have attracted considerable interest over the last few years. Besides 
than in four-dimensional models H], they naturally arise in the context of compactified 
Kaluza- Klein theories ifloTl as well as in superstring theories fllll . Moreover, it has been 
argued that an ALP with mass m ~ 10 -33 eV is a good candidate for quintessential 



dark energy [12] which presumably triggers the present accelerated expansion of the 
Universe. 

Specifically, ALPs are spin-zero light bosons defined by the following low-energy 
effective lagrangian 

^alp = ^ad^a-^m 2 a 2 -^F^F^a, (3) 

where F^ v is the electromagnetic field strength, Fny is its dual and a denotes the ALP 
field. According to the above view, it is assumed that the inverse two-photon coupling M 

— 1/2 

is much larger than G F ~ 250 GeV. On the other hand, the ALP mass m is supposed 

-1/2 

to be much smaller than G F ~ 250 GeV and for definiteness we take m < 1 eV. As far 
as generic ALPs are concerned, the parameters M and m are regarded as independent. 



The situation is at variance with the case of the standard Axion 111311 . for which the 
relationship 

10 10 GeV\ 



m = Q.lk-\ — leV (4) 



holds, with k a model-dependent constant of order one [|14fl 



General considerations 



A characteristic feature of ALPs is the trilinear y-y-a vertex described by the last 
term in JzfALP* whereby one ALP couples to two photons. This interaction gives rise 



to photon- ALP conversion, which leads in turn to a very interesting variety of physical 
processes, including the one to which the present paper is devoted. 

One of these processes consists in ALP photo-production through the Primakoff 
process, which takes place when an incoming photon scatters on a charged particle 
thereby becoming an ALP upon the exchange of a virtual photon. Hot, dense plasmas 
in stellar cores are ideal environments where the Primakoff process involving thermal 
photons can occur. Once produced, the ALPs escape because their mean free path is 
much larger than the stellar radius, thereby carrying off energy. Owing to the virial 
equilibrium, the core has a negative specific heat. Therefore it reacts to such an energy 
loss by getting hotter. As a result, the rate of nuclear reactions sharply increases, bringing 
about a substantial change in the observed properties of stars. Still, current models of 
stellar evolution are in fairly good agreement with observations. Hence, M has to be 
large enough to provide a sufficient suppression of unwanted ALP effects. This argument 
has been applied in a quantitative fashion to the Sun, to main-sequence stars and to red- 
giants stars, with the result jl5\ \ 

M>10 10 GeV. (5) 

Remarkably enough, for m < 0.02 eV almost the same conclusion, namely 

M>0.86-10 10 GeV, (6) 

has been reached by the negative result of the CAST experiment designed to detect ALPs 
emitted by the Sun Rial . 

Another process implied by the y-y-a vertex is direct photon- ALP conversion without 
the need of charged particles playing the role of catalysts, provided an external magnetic 
(electric) field is present. Whenever the external field extends over a large region and 
the momentum transfer is small, photon-ALP conversion becomes coherent and can 
be regarded as an oscillation phemomenon. Indeed, the y-y-a vertex produces an off- 
diagonal element in the mass matrix for the photon-ALP system in the presence of 
an external field. Therefore, the interaction eigenstates differ from the propagation 
eigenstates and photon-ALP oscillations show up lfl7ll . The situation is analogous to 
what happens in the case of massive neutrinos with different flavours, apart from an 
important difference. All neutrinos have equal spin, and so neutrino oscillations can 
freely occur. Instead, ALPs are supposed to have spin zero whereas the photon has 
spin one, and so one of them can transform into the other only if the spin mismatch 
is compensated for by an external field. 

One consequence of photon-ALP oscillations is that a lower bound on M stronger than 
condition © can be derived for ALPs with mass m < 10~ 10 eV. In this connections, 
two methods have been put forward. One concerns the energetics of the supernova 
1987a. Because the emitted neutrinos have been observed and the whole energy budget 
is rather well known, an upper bound on the ALP flux can be derived UM- Alternatively, 
observation of time-lag between opposite-polarization modes in pulsar radio emission 
similarly yields an upper bound on the two-photon coupling of an ALP |fl^| . In either 
case, one gets 

M>3-10 n GeV. (7) 



Application to VHE gamma-ray observations 



Our proposal now starts to become clear. We imagine that photons emitted by a distant 
blazar can oscillate into ALPs in the presence of cosmic magnetic fields. So, the emitted 
flux gets reduced along the line-of-sight because some photons become ALPs. If this 
were the whole story, the observable prediction would merely be a dimming Q20Q . Things 
can be totally different when photon absorption becomes important. To see how this 
comes about, suppose that a sizeable fraction of the emitted photons convert into ALPs 
close enough to the source and that a nonnegligible fraction of the ALPs in question are 
in turn converted back into photons close enough to the Earth. Because ALPs propagate 
unimpeded, in such a situation the observed photon flux can be considerably larger 
than the one predicted by Eq. ©. Below, we will turn this qualitative picture into a 
quantitative estimate. 

In principle, we have to evaluate the probability Py_>y(D) that a photon remains a 
photon after propagation from the source to us when allowance is made for photon- ALP 
oscillation as well as for photon absorption. 

However, an exact treatment would be exceedingly difficult, mainly because of the 
uncertainty concerning the configuration of the magnetic field responsible for photon- 
ALP oscillations. Actually, the line-of-sight to a distant blazar is expected to cross 
magnetic fields on different scales. A magnetic field is likely present in the source itself. 
In addition, the Galactic magnetic field can give a nontrivial contribution to the effect 
under consideration. Finally, extragalactic magnetic fields can play an important role. 
Throughout, we restrict our attention to extragalactic magnetic fields, whose existence 
has been demonstrated very recently by AUGER observations ||21|1 . A complementary 
picture involving only magnetic fields in the source and in the Milky Way has been 
considered in ref. Q22(] . 

Unfortunately, almost nothing is known about the morphology of extragalactic mag- 
netic fields, which reflects both their cosmic origin and the evolutionary history of bary- 
onic matter. While it is evident that their coherence length cannot be arbitrarily large, no 
reliable estimate of its value is presently available. As far as our analysis is concerned, 
this means that we cannot suppose that extragalactic magnetic fields are constants over 
the whole distance to the source - yet, their spatial dependence is unknown. The usual 
way out of this difficulty amounts to suppose that extragalactic magnetic fields B have 
a domain-like structure Q23[|. That is, B is assumed to be constant over a domain of size 
Ldom equal to its coherence length, with B randomly changing its direction from one 
domain to another but keeping approximately the same strength. Reference values to 
be used throughout are B ~ 10 -9 G and Ldom — 1 Mpc, which are close to existing up- 
per limits but consistent with them ll24ll . Remarkably enough, such a picture of cosmic 
magnetic fields turns out to be consistent with AUGER observations Q25Q. 

Thus, the whole propagation process of the photon beam can be recovered by iterating 
the propagation over a single domain as many times as the number of domains crossed 
by the beam, taking each time a random value for the angle 6 between B and a fixed 
fiducial direction. In this way, we are effectively led to the much easier problem of 
photon- ALP oscillation in a constant magnetic field. 

Another simplification is brought about by the fact that photon absorption is insen- 
sitive to the specific properties of the photon-ALP oscillation - it only depends on the 



photon mean free path Xy. So, we can first identify the regime in which Py_>y is maximal 
over a single magnetic domain and next we can work out how much this probability is 
suppressed by photon absorption. 

Consider in the first place the propagation of a photon beam over one domain without 
photon absorption. Since now B is constant, the probability that a photon converts to an 
ALP after a distance y can be computed exactly and reads |Q/z|] 

P^ a (y) = sin 2 2a sin 2 , (8) 

where the photon- ALP mixing angle a is 

1 / B \ 

a = - arcsin — — — (9) 
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and the oscillation wave number reads 
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so that the oscillation length is L osc = 27r/A osc . Actually, eq. (PTOl) pertains to the situation 
in which the beam propagates in a cold intergalactic plasma with plasma frequency 

«pi= -3.69-10- 11 ^3 eV, (11) 

V f^ig J > cm / 

where n e denotes the electron density. Because we are dealing with weak magnetic fields, 
their contribution to the vacuum refractive index is totally negligible. 

A look back at Eqs. d8]), © and (flOl) shows that the photon- ALP transition probability 
is maximal in the strong-mixing regime, namely for E \m 2 — ai^M/lB. A conserva- 
tive estimate of the density of intergalactic plasma yields n e ~ 10~ 7 cm~ 3 |26|], result- 
ing in the plasma frequency (Op\ ~ 1.17 • 10~ 14 eV. Therefore, the strong-mixing condi- 
tion takes the explicit form |(m/10" 10 eV) 2 - 1.37 • 10 8 | < 0.38(£/GeV)(fl/10~ 9 G) 
(10 10 GeV/M). Since we are interested in the energy range E > 10 2 GeV, we find that 
for e.g. M > 4 • 10 11 GeV the bound m < 10~ 10 eV has to be satisfied. Recalling Eq. ©, 
we see that the standard axion is excluded. We also remark that the present mechanism 
works for arbitrarily small values of m, provided M happens to be considerably smaller 
than the Planck mass Mp ~ 1.22 • 10 19 GeV. As a consequence, our result also applies 
to models of dark energy wherein quintessence enjoys a photon coupling described by 



-^alp II 1211 . thus ultimately providing an important observational test for these models. 



ENERGY SPECTRUM FOR BLAZAR 3C 279 



It is straightforward to solve the beam propagation equation once photon absorption by 
the EBL is taken into account and produces a finite photon mean free path X r [|23j]. In the 



strong-mixing regime, the probability for a photon to become an ALP after a distance 
y < Atom becomes 
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whereas the probability that a photon remains a photon now reads 
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where we have introduced the dimensionless parameter 

BXy 5 \ /10 n GeV\ / A* 
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The relation between the photon energy E and the source redshift z yields for 3C 279 
A y ~ 450 Mpc at E = 500 GeV so that we have 5 ~ 12.4 in this case study. 

Over distances y ^> L& om , the transition probabilities Py-^ a (y) and Py-^y(y) arise as 
the incoherent average of Eqs. (PT21) and (fT3l over /V ~ Cy/Atom) domains crossed by the 
beam, respectively. Assuming (as before) that the beam propagates along the y direction 
and choosing the x and z directions arbitrarily in the orthogonal plane, the problem 
becomes truly three-dimensional, because of the random orientation of the magnetic 
field. Consequently, the beam state is described by the vector (y x , y z ,a). 

We have derived the propagation equations describing the absorption of photons due 
to the interaction with the EBL and their oscillations to ALPs (and vice-versa). As in 
Ref. [23], we are led to the transfer equation 
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Here 6 denotes the angle between the x axis and the extragalactic B in a single domain. 
Starting with an unpolarized photon beam, we propagate it by iterating Eq. (fT5l) as many 
times as the number of domains crossed by the beam, taking each time a random value 
for the angle e (this reflects the random orientation of B). We repeat such a procedure 
10 000 times and finally we average over all these realizations of the propagation pro- 
cess. The resulting spectrum is exhibited in Fig.[2l We find that about 13% of the photons 
arrive to the Earth for E = 500 GeV, resulting in an enhancement by a factor of about 
20 with respect to the flux expected in the absence of the proposed oscillation mech- 
anism; the comparison is made with the "best-fit model" described in Kneiske et al. 
(2004) [2]. The same calculation gives a fraction of 76% for E = 100 GeV (to be com- 
pared to 67% without the oscillation mechanism) and a fraction of 3.4% for E = 1 TeV 
(to be compared to 0.0045% without the oscillation mechanism). In addition, we have 
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FIGURE 2. The two lowest lines give the fraction of photons surviving from a source at the same 
distance of 3C 279 without the oscillation mechanism, for the "best-fit model" of EBL (dashed line) and 
for the minimum EBL density compatible with cosmology [2]. The solid line represents the prediction 
of the oscillation mechanism for B ~ 10~ 9 G and Ld om — 1 Mpc within the "best-fit model" of EBL. The 
gray band is the envelope of the results obtained by independently changing B and Ldom within a factor of 
10 about their preferred values. 



checked the stability of our result against independent variations of B and Ld om within a 
factor of 10 about their preferred values. The corresponding spectrum is represented by 
the gray band in Fig. [21 We remark that the standard deviation of the above averaging 
procedure lies well inside the gray band. Our prediction can be tested with the satellite- 
borne Fermi/LAT detector as well as with the ground-based IACTs H.E.S.S., MAGIC, 
CANGAROO III, VERITAS and the Extensive Air Shower arrays ARGO-YBJ and MI- 
LAGRO. Our result also offers an important observational test for models of dark energy 
wherein quintessence is coupled to the photon through an effective dimension-five op- 
erator. 
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